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Abstract 

Continuing work initiated in an earlier publication [Yamada, Asada, Phys. Rev. D 82, 104019 
(2010)], we investigate collinear solutions to the general relativistic three-body problem. We prove 
the uniqueness of the configuration for given system parameters (the masses and the end-to-end 
length). First, we show that the equation determining the distance ratio among the three masses, 
which has been obtained as a seventh-order polynomial in the previous paper, has at most three 
positive roots, which apparently provide three cases of the distance ratio. It is found, however, 
that, even for such cases, there exists one physically reasonable root and only one, because the 
remaining two positive roots do not satisfy the slow motion assumption in the post-Newtonian 
approximation and are thus discarded. This means that, especially for the restricted three-body 
problem, exactly three positions of a third body are true even at the post-Newtonian order. They 
are relativistic counterparts of the Newtonian Lagrange points Li, L2 and L3. We show also 
that, for the same masses and full length, the angular velocity of the post-Newtonian collinear 
configuration is smaller than that for the Newtonian case. Provided that the masses and angular 
rate are fixed, the relativistic end-to-end length is shorter than the Newtonian one. 

PACS numbers: 04.25.Nx, 95.10.Ce, 95.30.Sf, 45.50.Pk 
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I. INTRODUCTION 



The three-body problem in Newtonian gravity belongs among classical problems in as- 
tronomy and physics (e.g, [l, 2|). In 1765, Euler found a coUinear solution for the restricted 
three-body problem, where one of three bodies is a test mass. Soon after, his solution was 
extended for a general three-body problem by Lagrange, who also found an equilateral tri- 
angle solution in 1772. Now, the solutions for the restricted three-body problem are called 
Lagrange points Li, L2, L3, L4 and L5, which are described in textbooks of classical mechan- 
ics SOHO (Solar and Heliospheric Observatory) and WMAP (Wilkinson Microwave 
Anisotropy Probe) launched by NASA are in operation at the Sun-Earth Li and L2, re- 
spectively. LISA (Laser Interferometer Space Antenna) pathfinder is planned to go to Li. 
Lagrange points have recently attracted renewed interests for relativistic astrophysics 
where they have discussed the gravitational radiation reaction on L4 and L5 analytically |4| 

and by numerical methods aaa. 

As a pioneering work, Nordtvedt pointed out that the location of the triangular points 
is very sensitive to the ratio of the gravitational mass to the inertial one jl]. Along this 
course, it is interesting as a gravity experiment to discuss the three-body coupling terms 
at the post-Newtonian order, because some of the terms are proportional to a product of 
three masses as Mi x M2 x M3. Such a term appears only for relativistic three (or more) 
body systems: For a relativistic binary with two masses Mi and M2, there exist MfM2 and 
MiM| without such a three-mass product. For a Newtonian three-body system, we have 
only the two-body coupling terms proportional to M1M2, M2M3 or M3M1. 

The relativistic perihelion advance of Mercury is detected only after much larger shifts due 
to Newtonian perturbations by other planets such as the Venus and Jupiter are taken into 
account in the astrometric data analysis. In this sense, effects by the three body coupling 
are worthy to investigate. Nevertheless, most of post-Newtonian works have focused on 
either compact binaries because of our interest in gravitational waves astronomy or N-body 
equation of motion (and coordinate systems) in the weak field such as the solar system 

n hn 

(e.g. [8]). Actually, future space astrometric missions such as Gaia [9|, llO[ require a general 
relativistic modeling of the solar system within the accuracy of a micro arc-second 11 |. 



Furthermore, a binary plus a third body have 
gravitational waves induced by the third body 



)een discussed also for perturbations of 



12 



15|. 
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After efforts to find a general solution, Poincare proved that it is impossible to describe 
all the solutions to the three-body problem even for the 1/r potential. Namely, we cannot 
analytically obtain all the solutions. Nevertheless, the number of new solutions is increasing 



16| . Therefore, the three-body problem still remains an open issue even for Newton gravity. 



The theory of general relativity is currently the most successful gravitational theory 
describing the nature of space and time. Hence, it is important to take account of general 
relativistic effects on three-body configurations. The figure-eight configuration that was 
found decades ago [l7, 1^ has been recently studied at the first post-Newtonian [lol and 
also the second post-Newtonian orders [20] . According to their numerical investigations, the 
solution remains true with a slight change in the figure-eight shape because of relativistic 
effects. 

On the other hand, the post-Newtonian coUinear configuration obtained in the previous 
paper [2l| may offer a useful toy model for relativistic three-body interactions, because it is 



tractable by hand without numerical simulations. This solution is a relativistic extension of 
Euler's coUinear one, where three bodies move around the common center of mass with the 
same orbital period and always line up. 



In fact, their formu 
ratio among masses 



ation leads to a seventh-order equation determining the distance 
2l| . Here, it should be noted that only positive roots are acceptable, 
because the distance ratio must be positive. Properties of the master equation have not been 
known yet. How many positive roots for it are there? The main purpose of this paper is to 
analytically investigate the number of the positive roots. In particular, we shall prove the 
uniqueness of the configuration for given system parameters (the masses and the end-to-end 
length) . 

This paper is organized as follows. In section II, we briefly summarize formulations for 
coUinear solutions at the Newtonian and post-Newtonian orders. We discuss positive roots 
for the seventh-order equation for determining the distance ratio in section HI. In section IV, 
we show the uniqueness of the configuration for given system parameters (the masses and the 
end-to-end length). We also compare the angular velocity of the post-Newtonian coUinear 
configuration with that for the Newtonian one. Section V is devoted to the conclusion. We 
provide some detailed calculations regarding the angular velocity of coUinear configurations 
in the Appendix. 

Throughout this paper, we take the units of G = c = 1. 
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II. EQUATION FOR THE DISTANCE RATIO AMONG THREE MASSES 



Let us begin by summarizing the derivation of the Euler's colhnear solution for the circular 
three-body problem in Newton gravity. We consider Euler's solution, for which each mass 
moves around their common center of mass denoted as Xq with a constant angular velocity 
LO. Hence, it is convenient to use the corotating frame with the same angular velocity uj. We 
choose an orbital plane normal to the total angular momentum as the x — y plane in such 
a corotating frame. We locate all the three bodies on a single line, along which we take the 
x-coordinate. The location of each mass M/ (/ = 1, 2, 3) is written as Xi = {xj, 0). Without 
loss of generality, we assume x^, < X2 < Xi. Let Ri define the relative position of each mass 
with respective to the center of mass Xq = {xq, 0), namely Rj = xj —xq {Ri |^/| unless 
a^G = 0). We choose x = Q between Mi and M3. We thus have R3 < R2 < Ri, R3 < and 
Ri > 0. 

It is convenient to define a ratio as R23/R12 = z, which is an important variable in the 
following formulation. Then we have R13 = (1 + z)Ri2- The equation of motion becomes 

2 M2 M3 
= -^ + ^' (1) 

-"-12 ^IZ 

^ Ml My, 



p 2 Ml M2 

-"-13 -'^23 



where we define 



Rij = Xi- Xj, (4) 
Rij = \Rij\. (5) 

First, we subtract Eq. ([2]) from Eq. ([1]) and Eq. ([3]) from Eq. ([2]) and use -R12 = |Xi — X2I 
and R23 = 1X2— X^l. Such a subtraction procedure will be useful also at the post-Newtonian 



order, because we can avoid directly using the post-Newtonian center of mass [23 



2^. Next, 



we compute a ratio between them to delete cu^. Hence a fifth-order equation is obtained as 

(Mi+M2)z' + (3Mi+2M2)/ + (3Mi+M2)z='-(M2+3M3)2'-(2M2+3M3)2-(M2+M3) = 0. 

(6) 



Now we have a condition as 2; > 0. Descartes' rule of signs (e.g., 22]) states that the number 



of positive roots either equals that of sign changes in coefficients of a polynomial or less than 



it by a multiple of two. According to this rule, Eq. ([6]) has only the positive root z > 0, 
though such a fifth-order equation cannot be solved in algebraic manners as shown by Galois 



(e.g. 



22|). After obtaining z, one can substitute it into a difference, for instance between 
Eqs. ([T]) and (E]). Hence we get u. 

In order to include the dominant part of general relativistic effects, we take account of 
the terms at the first post-Newtonian order. Namely, the massive bodies obey the Einstein- 



Infeld-Hoffman (EIH) equation of motion as 23l. \24\ 



dv 



K 



dt 
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Ma 
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RakR^ 



(7) 



CA 



A+K Cj^A 

where Vi denotes the velocity of each mass in an inertial frame and we define 

Rij 



nij = 



R 



ij 



(8) 



and we assume the slow motion {\vj\ ^ c). 

We obtain a lengthy form of the equation of motion for each body. By subtracting the 
post-Newtonian equation of motion for M3 from that for Mi for instance, we obtain the 
equation as jsil 

Ri3U^ = Fn + Fm + Fvu\ (9) 

where we denote a = -R13 and the Newtonian term and the post-Newtonian parts Fm 
(dependent on the masses only) and Fy (velocity-dependent part divided by w^) are defined 
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as 



M 



[U, + U^)z' + (1 - Z/i - Z/3)(l + + zf 

{A-Aui + z/3)(l - z/1 - z/3) 
+ (12-7z/i + 3z/3)(l-z/i-z/3);2 

+ (12-Z/i + Z/3)(l-Z/i-Z/3)z2 

+ (8 - 7z/i - 7z/3 + 8z/ii/3 + 3i^i + 31/1)^3 

+ (12 + 1/1 -Z/3)(l-I/i -Z/3)Z^ 

+ (12 + 3z/i-7z/3)(l-z/i-z/3)2^ 

+ (4 + Z/i-4z/3)(l-Z/i-Z/3)/ 



(10) 



(11) 



M 



1 + zfz^ 



-ulil - z/1 - Z/3) 

— 2z/i(l + z/1 — z/3)(l — z/1 — z/s)^; 

+ (2 - 2z/i + z/3 + 6Z/1Z/3 - 3Z/3 + z/j' - 3z/^z/3 - 3Z/1Z/3 + z/g)^^ 

+2(2 - z/i - z/3)(l + z/i + z/3 - z/^^ + Z/1Z/3 - z/|)2;^ 

+ (2 + z/i - 2z/3 - 3z/i^ + 6Z/1Z/3 + vl - 3z/^z/3 - 3z/iz/| + z/|)2^ 

-2Z/3(1 - z/1 + Z/3)(l - z/1 - Z/3)25 

-z/|(l - z/1 - v^)z^ 



(12) 



respectively. Here, we define the mass ratio as z// = Mi/M for tlie total mass M = Yli^i 
and make a frequent use of 1/2 = 1 — i^i — ly^. It should be noted that in this truncated 
calculation we ignore the second post-Newtonian (or higher order) contributions so that we 
can replace, for instance, fi by Riu (using the Newtonian Ri) in post-Newtonian velocity- 
dependent terms such as vf. 

In a similar manner to the above Newtonian formulation, straightforward but lengthy 
calculations lead to a seventh-order equation as |2l[ | 

7 

F{z) = J2^kz'' = 0, (13) 

fe=0 
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where we define 

A7 = 



M 
a 



-4 - 2(z/i - 4Z/3) + 2{ul + 2Z/1Z/3 - 2z/|) - 2uiU3{ui + z/3) 
M 



1-Z/3 + 



-13 - (lOz/i - I7Z/3) + 2(2//!^ + 8z/iz/3 - u. 



= 2 + 1^1- 2u-i 



M 



Ai = 1 + 2z/i - z/3 + 



M 



-15 - (18z/i - 5z/3) + 4(5z/iZ/3 + 4z/| 



+6(z/f - z/iz/f - ul) 



-6 - 2(5z/i + 2Z/3) - 4(2//!^ - Z/1Z/3 - 4z/|) 



+2(3z/f + u^us - 2vi1jI - 3z/| 



A3 = -(1 - z/i + 21/3) + 



M 



A, 



_(2 - 2z/i + z/3) + 



M 



6 + 2(2z/i + 5z/3) - 4(4z/2 + z/iZ/3 - 2z/| 



+2(3z/i^ + 2z/i2z/3 - z/iz/| - 3z/i^ 



15 - (5z/i - I8Z/3) - 4(4z/2 + 5Z/1Z/3) 



+6(z/i3 + z/2z/3 - vl) 



A = -(l-Z/i) + 



M 



13 - (17z/i - lOz/3) + 2(z/2 - 8z/iZ/3 - 2z/|) 



+2(z/f + 3z/2z/3 + 2uiul - ul] 



An = — 



4 - 2(4z/i - z/3) + 2(2z/2 - 2Z/1Z/3 - z/|) + 2z/iz/3(z/i + z/3) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



Here, the sign of Eq. ( !2T]) is chosen so that it can agree with the fifth-order equation Eq. 
(E]) in the Newtonian limit of M/a — > 0. This seventh-order equation is antisymmetric for 
exchanges between z/i and z/3, only if one makes a change as z ^ 1/z. This antisjTiimetry 
may vahdate the complicated form of each coefficient. 
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Once a positive root for Eq. (fT3|) is found, the root z can be substituted into Eq. ([9]) in 
order to obtain the angular velocity cu. 

The angular velocity including the post-Newtonian effects is obtained from Eq. Q as 



where a; at = (F/v/i^is)^/^ denotes the angular velocity of the Newtonian collinear orbit. Note 
that the slow motion is assumed to derive Eq. (1221) which is analogous to Kepler's third law. 

III. EXISTENCE OF POSITIVE ROOTS 

In this section, we show that there always exist positive roots for the seventh-order 
equation that has been derived as Eq. f|T3|l . This is nothing but the existence of the post- 
Newtonian collinear solution. 

For later convenience, we recover and thus rewrite a coefficient Aq as 



which immediately leads to > 0. 

In a similar manner, one can show < 0. An alternative but powerful way to see this is 
using the antisymmetry of the seventh-order equation for transformations between masses 
Ml and M3 as v\ o 2/3 and 2; ^ This transformation makes a change as — )■ —A^. 
By using A^ > 0, therefore, we have always Aj < 0. 

Bringing the above results together, we have F{0) = Aq > and F(oo) = Ajz'^\z^oo < 0. 
Therefore, the number of positive roots for F{z) = either equals to one or more than it 
by a multiple of two. 

Let us investigate the seventh-order equation in order to more precisely determine the 
number of positive roots. We decompose each coefficient Aj. into the Newtonian part Aj^ik 
and the post-Newtonian one Apj^i^. Note that Ajvfc agrees with the coefficient of 2;^^^ (but 
not z'^) in Eq. ([6]). In the Newtonian fifth-order equation by Eq. ([6]), we have Anq > 0, 
An5 > 0, An4, > 0, Ans < 0, An2 < 0, Ani < 0. In the post-Newtonian approximation, the 
post-Newtonian parts must be much smaller than the Newtonian ones (|v4pArfc| ^ \ANk\ for 
each k), so that the post-Newtonian correction cannot change the sign of each coefficient A^. 
We thus have Aq > 0, ^5 > 0, A4 > 0, ^3 < 0, A2 < 0, Ai < 0. By combining them with 




(22) 



Aq = 2 (z/2 + Z^3)(2z^2 + 2^3 + ^2^3), 



(23) 
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TABLE I: Values of z, lo and aui for Figure [H Here are three positive roots, where we assume 

Zl < Z2 < Z3. 







Z2 


^3 


z 


3.635x10"'^ 


1.000 


2751 




8.723x10"^ 


2.449x10-*^ 


8.723x10"^ 


auj 


0.8723 


0.02449 


0.8723 



^7 < and y4o > 0, the number of sign changes of the coefficients in Eq. f|T3|l is necessarily 
three. Therefore, Descartes' rule of signs indicates that Eq. (fT3|) has either one or three 
roots. We can easily understand that one of them is a correction to the Newtonian orbit. 
What are the other two roots? We shall investigate them in next section. 

IV. UNIQUENESS OF THE POST-NEWTONIAN COLLINEAR SOLUTION 

Figure [T] shows that the equation has three positive roots, where we assume vi = 1/7, 
1^2 = 5/7, z/3 = 1/7, a/M = 10^ [v ~ 10~^). Table [T] shows numerical values of z, u and 
au for Figure [H Two out of the three roots do not satisfy a slow-motion condition for the 
post-Newtonian approximation as shown below. 

Here we show that the remaining two positive roots must be discarded. Because of the 
antisymmetry of Eq. f[T^ for the transformation as z ^ I/2;, the two roots must be a pair 
through this transformation associated with exchanges between Mi and M3. Let the smaller 
root and the larger one be denoted as zs and zl, respectively. 

First, we consider the smallest positive root zs, where we assume zs '^l- Then, Eq. f|T3|) 
is approximated as 

Aizs + Ao = 0, (24) 

where Aq starts at the post-Newtonian order without Newtonian terms and Ai = Aj^i+Apj^i 
has both the Newtonian terms and post-Newtonian corrections (I^Aril ^ |y4pAri|). We thus 
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obtain an approximate form of the smallest root as 

^0 



Am 

O {'A . (25) 



\ a 

where we used Eqs. ( l20i) and ( l2Ti) . This implies that zs is indeed of the post-Newtonian 
order, in consistent with zs <^ 1. At this point, however, we cannot discard this smallest 
root Zs- 

As a next step, let us make an order-of-magnitude estimation for the angular velocity ujs 
that satisfies Eq. (|9]) for zs, where us denotes the angular velocity corresponding to zs- We 
obtain from Eqs. (fTO|). (ITT]) and (fT2|) 



o 



M 



0(^1, (26) 



o 



M 
M' 



a^z\ 



o 



M 
If 



We thus find -R13 = a <C -Fy because M ^ a. Therefore, we find ~ Fm ~ Fyuj\ ^ -Ris^l 
in Eq. ([9]). This leads to 

ujs = O , (29) 

though bj\; = 0{M/a^) for the Newtonian case. Eq. (129|) implies an extremely fast rotation, 
since the rotational velocity becomes vs ~ acos = 0(1), namely, comparable to the speed 
of light. This unacceptable branch of such an extremely fast motion contradicts the post- 
Newtonian approximation and does not satisfy Eq. fl2^ . Hence, zs must be abandoned. 

Next, we consider the largest positive root zl, where we assume zl ^ 1. Then, Eq.(fT3ll 
is approximated as 

AjZl + Ae = 0, (30) 
10 



where starts at the post-Newtonian order without Newtonian terms and Aq = A^e+ApN^ 
has both the Newtonian terms and post-Newtonian corrections (IAatqI ^ I^PTveD- We thus 
obtain an approximate form of the largest root as 

A7 

where we used Eqs. (IT^ and (ITS!) . This imphes that 2;^^ is indeed of the post-Newtonian 
order, in consistent with zl ^ 1. At this point, however, we cannot discard this largest root 

As a next step, let us make an order-of-magnitude estimation for the angular velocity ul 
that satisfies Eq. for zi, where ul denotes the angular velocity corresponding to zl- We 
obtain from Eqs. ([TU]). (ITTj) and ([T^ 



o 



a2 



o 



M^zj 



(33) 



Fv = O (Mzl) 

We thus find i?i3 = a <^ Fy because M ^ a. Therefore, we find F/v ~ -Fm ~ Fv^t ^ Ri3<^l 
in Eq. (|9]). This leads to 

= O (i) . (35) 

though = 0{M/a^) for the Newtonian case. Eq. (!35|) implies an extremely fast rotation, 
since the rotational velocity becomes vl ~ owl = 0(1), namely, comparable to the speed 
of light. This unacceptable branch of such an extremely fast motion contradicts the post- 
Newtonian approximation and does not satisfy Eq. ( l22l) . Hence, also zl must be abandoned. 
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We should remember the transformation as z •<-> I/2;, namely 1 •<-> 3. Hence, zg and 
zi, correspond to each other as zg — '^/zl- In this sense, it seems natural that the above 
argument for discarding z^ is very similar to that of zs- 

As a result, two of the three positive roots are discarded as unphysical ones. Hence, we 
complete the proof of the uniqueness. 

We mention an application of the uniqueness theorem for the restricted three-body prob- 
lem. We have three possibilities for choosing a test mass as Mi =0, M2 = or M3 = 0. For 
each case, we have only the single coUinear solution. Therefore, the three equilibrium points 
exist along the symmetry axis of the system, and they are a generalization of Lagrange 
points Li, L2 and L3. 

Before closing this section, we mention an interesting property of the angular velocity 
of the coUinear configurations. For the same masses and full length, we have always an 
inequality as 

uj < un, (36) 

which means that the post-Newtonian orbital period measured in the coordinate time is 
longer than the Newtonian one. Provided that the masses and angular rate are fixed, the 
relativistic length a is shorter than the Newtonian one. Detailed calculations are given in 
the Appendix. 

V. CONCLUSION 

Wc proved the uniqueness of the coUinear configuration for given system parameters (the 
masses and the end-to-end length). It was shown that the equation determining the distance 
ratio among the three masses, which has been obtained as a seventh-order polynomial in 
the previous paper, has at most three positive roots, which apparently provide three cases 
of the distance ratio. It was found, however, that there exists one physically acceptable root 
and only one. The remaining two positive roots are discarded in the sense that they do not 
satisfy the slow motion ansatz in the post-Newtonian approximation. 

Especially for the restricted three-body problem, exactly three positions of a third body 
are true even at the post-Newtonian order. They are relativistic counterparts of the New- 
tonian Lagrange points Li, L2 and L3. 

It was shown also that, for the same masses and full length, the angular velocity of 
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the post-Newtonian collinear configuration is smaller than that for the Newtonian case. 
Provided that the masses and angular rate are fixed, the relativistic length a is shorter than 
the Newtonian one. 

Our way of discussion seems to work at the second (and higher) post-Newtonian orders, 
because the slow motion approximation is a key in the above proof. Therefore, the unique- 
ness of collinear configurations for a three-body system may be true even at higher orders, 
precisely speaking, if the configuration has the Newtonian limit. It is an open question 
whether fully general relativistic systems admit a particular solution that can appear only 
for a fast motion case and thus has no Newtonian limit. 

We would like to thank the referee for useful comments on the earlier version of the 
manuscript. We are grateful to Y. Kojima for useful conversations. This work was supported 
in part (H.A.) by a Japanese Grant-in-Aid for Scientific Research from the Ministry of 
Education, No. 21540252. 

Appendix A: Detailed calculations on the angular velocity 

Let us prove a; < c^Af- Eq- f l22|) is rewritten as 



Here, Fn is positive and thus the denominator of the R.H.S. of Eq. f lAip is always positive. 
What we have to do is to investigate the sign of the numerator for the R.H.S. of Eq. (lAip . 
The numerator is factored as 



u — un _ FmRiz + FvFn 



(Al) 




(A2) 
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where we define 

oio = -(1 - 1^1 - i^sfi^i, (A3) 
og = -(1 - i/1 - i/3)(4 + i/i - 2i/3 - 4i/ii/3 + 2i/^ + 21/^^3 - 2i/ii/| - 4z/|), (A4) 

as = -(1 - i^i - i^3)(18 + 4i/i - 9i/3 + 3i/i - 14i/ii/3 

+2i/|-z/j' + 7z/iZ/3 + 2z/iZ/|-6z/|), (A5) 

07 = -(1 - i/i - 1/3) (32 + 4i/i - 13i/3 + 12ul - 18i/ii/3 

+10i/| - Apf + 81/1^1/3 + 4i/ii/| - 8z/|) , (A6) 

ag = -(30 - 30ui - 37i/3 + Idu^ - 12uiUs + 27i/| - 22uf + 18^1^3 

+12i/ii/| - 28i/| + 6i/f - 4i/fi/3 - 15i/ii/| + 6i/ii/| + (A7) 

05 = -2(12 - 13i/i - 13i/3 + lli/i - lOi/ii/3 + lli/| - lli/f + 17iyfu3 

+17iyiui - lli^l + 4i/f - Sufus - l^vlvl - ^v^vl + 4i/^), (A8) 

04 = -(30 - 37i/i - 30i/3 + 27i/i - 12i/ii/3 + 19i/| - 28i/f + 12i/ii/3 

+I81/11/3 - 22vl + 111/^+ 6i/fi/3 - Ibvl^l - Av-^vl + Gi/g), (A9) 

03 = -(1 - i/i - i/3)(32 - 13i/i + 4i/3 + lOi/i - 18i/ii/3 

+12i/| - 8i/f + 4i/^i/3 + 8i/ii/| - 4i/|), (AlO) 

a2 = -(1 - i^i - i/3)(18 - Qi^i + 4i/3 + 2i/i - 14i/ii/3 

+3i/|-6i/f + 2i/^i/3 + 7i/ii/|-i/|), (All) 

ai = -(1 - i/i - i/3)(4 - 2i/i + i/3 + 2i/i - 4i/ii/3 - Avl - 2^1^^ + 2vivl), (A12) 



ao = -(l-i/i-i/3)V. 1, (A13) 



We show Ofc < for each k. It is trivial that qq < and aio < 0. For z/i -H- 1/3, we have 
a symmetry between ag ^ ai, ag ^ 02, 07 ^ 03 and ^ a^^. Therefore, it is sufficient to 
examine Og, ag, 07, and 05. 

First, let us show ag < 0. The nontrivial factor in Eq. flA4p turns out to be positive by 
noting the following relation as 

4 - 2z/3 - 4z/iz/3 + 2ul - 2z/iz/f - 4z/| 

4 2 
= 4(1 - v^f + 10z/3(l - — z/i - z/3 - Y^^'i^'s) 

6 

> 4(1 - u^f + 10z/3(l - —ui - Z/3) 

> 4(l-Z/3)3 + 10z/3(l-Z^l-/^3) 

> 0, (A14) 

where we used z/i ^ 0, 1/3 ^ 1, z/iZ/3 ^ z/i. Hence we find Og < 0. 

We discuss the sign of ag. By using ui = \ — U2 — to delete vi and recover z/2, the 
R.H.S. of Eq. (lASp is factored as 



- z/2 (24 - 7z/2 - 23z/3 + 4z/| + + lOz/^z/3 + 15z/2Z/|). (A15) 

One can show that the latter three terms are positive, since 

24 - 7z/2 - 23z/3 = 1 + 7(1 - z/2 - z/3) + 16(1 - z/3) 

> 0. (A16) 



Hence, the second factor in Eq. (lAlSP is always positive, which leads to ag < 0, and also 
02 < 0. 



Next, we examine 07. By recovering z/2 to delete z/i, the R.H.S. of Eq. (]A6ll is factored as 

- z/2(44 - 16z/2 - 39z/3 + 2Z/2Z/3 + 16z/| + 4z/^ + 2OZ/2Z/3 + 24z/2Z/|). (A17) 

A key thing is a positive as 

44 - 16z/2 - 39z/3 = 5 + 16(1 - z/g - y^) + 23(1 - z/3) 

> 0, (A18) 

which immediately leads to 07 < 0, and also 03 < 0. 
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We investigate qq. Similarly to ay, it is factored as 
- (3 + 34z/2 - z/3 - llz/2 - 34Z/2Z/3 + z/| - + 24z/2Z^| + 6z/2 + 28z/^z/3 + 331^2^^!). (A19) 
Here, we see that the following two combinations both are positive, 

34z/2 - llz/2 " 34z/2Z/3 = llz/2(l - z/2 - ^'s) + 23^2(1 - z/3) 

> 0, (A20) 
3-z/3-2z/^ > 2-2z/^ 

> 0, (A21) 



which show that Eq. ( ]A19I) is always negative. Hence, we show < 0, and also 04 < 0. 



Also for 05, it is factored as 



-2(3 + 8z/2 - z/3 + 2z/2 - IIZ/2Z/3 + z/| - 5z/| 

-7z/^z/3 + 12z/2Z/| + 4z/^ + 19z/^z/3 + 19z/^z/|). (A22) 

One can find the following rather tricky manipulation as 

3 + 8z/2 — z/3 + 2z/2 — llz/2^3 — 5z/2 — 7Z/2Z/3 

= 3 + 5Z/2[1 - Z/2(z/2 + z/3) - z/3] - z/3 - 6Z/2Z/3 

> 3 + 5z/2(l - z/2 - z/3) - z/3 - 6z/2Z/3 

> ^, (A23) 

where we used < z/2 + z/3 < 1 and Z/2Z/3 ^1/4. Hence, we find 05 < 0. 

As a consequence, all the coefficients are always negative, which shows oj < ujn for any 
mass ratio. 
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FIG. 1: Top panel: The seventh-order polynomial in the L.H.S. of Eq. (|13p . The horizontal axis is 
chosen as z. We take vx = 1/7, = 5/7, 1^3 = 1/7, a/ M = {v ^ IQ-^) in order to exaggerate 
small effects in these figures. Clearly such a symmetric choice of the mass ratios produces a trivial 
root as z = 1, which makes it easy to check numerical calculations. M2 is relatively large so that 
the centrifugal force can be large. 

Middle panel: The seventh-order polynomial around the smallest positive root zs- 
Bottom panel: The polynomial around the moderate positive root. 
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